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Abstract 

Based on the concept of complementary media, we propose a novel design which can enhance the 
electromagnetic wave scattering cross section of an object so that it looks like a scatterer bigger 



^ . than the scale of the device. Such a "superscatterer" is realized by coating a negative refractive 



material shell on a perfect electrical conductor cylinder. The scattering field is analytically obtained 
by Mie scattering theory, and confirmed by full-wave simulations numerically. Such a device can 
be regarded as a cylindrical concave mirror for all angles. 

PACS numbers: 41.20.Jb, 42.25.Fx 



Recently great progress has been 
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nipulating the electromagnetic (EM) fields by means of metamaterials. By employing the 
coordinate transformation approach proposed by Leonhardt [l| and Pendry et al. 121, various 

nnn 

exciting functional EM devices have been reported [^, |9|, HOj. This methodology provides 
a clear geometric picture of those designed devices, and the permittivity and permeability 
tensors of designed functional materials can be derived from coordinate transformations di- 
rectly. On the other hand, Mie scattering theory |lll . Il2l . Il4| provided an analytic approach 
to quantitatively analyze the scattering properties of EM fields, and the aforementioned 
functional devices can also be realized by choosing different scalar transformation functions. 
The combination of the geometric and the analytic approaches may give novel and surprising 
results. For example, in contrast to invisibility cloak, we may design an EM transformation 
media device to enlarge the scattering cross-section of a small object. In this way the object 
is effectively magnified to a size larger than the object plus the device so that it is much 
easer for EM wave detection, which we refer hereafter as a superscatterer. 

In this Letter, we propose a generalized technique to achieve such a superscatterer. In 
the quasistatic limit, Nicorovici et al. [l5| had demonstrated that the properties of a coated 
core can be extended beyond the shell into the matrix. It is called a partially-resonant sys- 
tem in which at least one dielectric constant of three-phase composite structure is negative. 



Moreover, by using the negative refracting material (NRM) ly, \\J\ , Pendry and Ramakr- 



ishna 181] showed how to image an object by perfect cylindrical lens based on the concept of 
complementary media. Here we limit our discussion on the 2D case and prove that a special 
magnified image acts as a real object for EM wave detection. 

Figure [1] shows a heuristic model stemmed from the concept of complementary media 
[l8|. In Fig. [T^a), a slab of metamaterial with e = /i = — 1 images point source 5* at F2 
inside and Fi outside the metamaterial if the distance between source and slab is less than 
the slab thickness. Now a perfect electrical conductor (PEC) boundary is set at -Ri, so the 
light propagating in the slab is completely refiected by it and focuses at I2 inside and Ji 
outside the slab. Since the slab and the vacuum with the same thickness are complementary 
to each other, it looks as if the medium between R^ and -Ri is moved, and the PEC at Ri is 
virtually shifted to R^ and takes effect as a real PEC boundary. The picture can be extended 
to the 2D case as shown in Fig. [U^b), in which z axis is perpendicular to the paper. In 
the 2D case, the i?i, R2 and -R3 represent the radius of the inner region, the outer radius of 





(b) 



FIG. 1: (a) The schematic demonstration of the behavior of a beam propagates in a complementary 
media and PEC boundary, (b) The behavior is extended to the 2D case, where a superscatterer is 
formed with its effective size shown by the dashed hne. 

the cylindrical annulus (metamaterial shell) and the effective radius of the virtual cylinder, 
respectively. In order to move effectively the inner PEC boundary at r = i?i to the surface 
of the virtual cylinder, where the radius is r = R^ > R2, both permittivity and permeability 
tensors in the metamaterial shell should be selected properly. The Mie scattering theory is 
a powerful tool to accomplish the goal. 

We consider a transverse-electric (TE) polarized EM incident field with harmonic time 
dependence exp {—iwt). The coordinate system is the cylindrical coordinate coinciding with 
the cylinder we considered. In this coordinate system the permittivity and permeability 
tensors can be put in the following general form 



er{r)ff + eg{r)99 + ez{r)zz, 



jj, = jjr{r)rr + iJ,Q{r)UU + jjz[r)zz. 



(1) 



The EM fields in the homogeneous material region (r > R2) are well known and can be 

expressed by the superposition of cylindrical functions. The wave equation of E^ in the shell 

is written as 

1 1 r) / T riP. \ ^ ^ r) /' ^ fiF. \ 

(2) 
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where fco is the wave vector of the EM wave in vacuum. Now we introduce a new coordinate 
system {f{r), 6, z), in which /(r) is a continuous and piecewise differentiable function of the 
original radial coordinate r, the wave equation of field Ez{f{r), 9) = Ez{r, 9) is transformed 



to 



''-f'4.i-f''-^]^^^ii^'4^]^^iE.^o, 



(3) 



where f'jr) denotes d/(r)/dr. By taking the components of permittivity and permeabihty 
tensors as (l4l |. 
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(4) 



where eo and /io are the vacuum permittivity and permeabihty, the Eq. ([3]) can be solved 
by separation of variables E^ = R{f)Q{6). The solution of R{f) and O(^) are just the 
r/ith-order Bessel functions and exp{im6), in which m are integers. 

With the above analysis, the electric fields in each domain are expressed as 



EJr^ 



0, 



r < Ri, 



Em (aUmikofir)) + al,H^rk\kof{r))) exp(^m^), R^ < r < R2, (5) 



Em {(^LJmikor) + (3^H}n'{kor) exp(imi 



(1)^ 



r > R2. 



r(i) 



where Jm and Hm are the mth-order Bessel function and Hankel function of the first kind, 
respectively. By means of the orthogonality of explimO) and the continuity of E^ and Hg 
at two interfaces (r = Ri and r = R2), we obtain the linear relationships between the 
coefficients as follows, 



a: 



a. 



,Jm{kof{Ri)) + a'mH^^\k,f{Ri)) = 0, 



(6a) 



R2 



f{R2) 



JUkof{R2)) + alH^^\kof{R2)) = PlnMkoR2) + P'mHi^\koR2), (6b) 
[alJiikofiR2)) + <ifi^)'(A;o/(i?2))] = (3inJ'mikoR2) + /3:,H^nkoR2), (6c) 



where the prime denotes differentiation with respect to the entire argument of Bessel func- 
tions. By imposing the boundary condition of /(r), /(-R2) = R2, we get the solutions of 
Eqs.(l6]) as follows. 






m = 0,±l,±2,.... 



(7) 



This is the exact solution for scattering matrix. 

We can draw some interesting conclusions from Eqs. ([7]). If /(r) is a monotonic function 
and f{Ri) = R3 > R2, the material of the shell must be NRM from Eqs.(jl]), and the fields 
at point {r,6,z) in inner annulus Ri < r < R2 are equal to those at point {f{r),6,z) in 
outer region R2 < r < R^. Moreover, in the region r > R^, the scattering fields are exactly 
equal to those scattered by a PEC cylinder with radius R3, as if whole cylindrical annulus 
of i?i < r < i?3 is moved and the PEC boundary at r = i?i is magnified and shifted to 
r = R3. Here, the cylinder with radius f{Ri) is called a virtual cylinder. It is just the 
expected result in Fig. [D^b). 




FIG. 2: A simple function /(r) satisfies the condition of FiglTJb) 



For simplicity, we choose a linear function for /(r) 



fir) 



bo{R2 - r)/(i?2 - Ri) + R2, Ri<r< R2, 
r, r > i?2- 



(8) 



Here 60 is an important parameter to determine the magnification factor. It should be 
noted that both cloak and concentrator can also be obtained in the same manner. For 
example, it is a perfect cloak when 69 = —R2, and it becomes an imperfect cloak when 
Ri — R2 > bo > —R2, or a concentrator when Bq > Ri — R2 and /(r) = rf{Rx)/Ri in the 
region r < Ri. However, the case 60 > 0, which provides a folded geometry [19|, l20|, is less 
discussed. In fact, with the help of the geometric picture of the concentrator, it means a 
big cylinder with radius 60 + R2 is compressed into a small cylinder with radius Ri , and the 
gap between r = Ri and r = Bq + R2 has to be filled by a pair of complementary media: 
One is the vacuum in R2 < r < bo + R2 and the other is the NRM shell in Ri < r < R2. 
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FIG. 3: Snapshot of the total and scattering electric field, (a) - (b) The total and scattering 
electric fields induced by PEC cylinder with radius R^ = O.Sm, respectively, (c) - (d) The total 
and scattering electric fields induced by the designed device (the radius of virtual cylinder is O.Sm), 
respectively. 

This device scatters the same fields as the uncompressed cylinder, which extends beyond 
the shell, so we call it "superscatterer" . 

Next, we present the patterns of electric field calculated by finite element solver of the 
Comsol Multiphysics software package. In Fig. [3l the plane wave is normal incident from left 
to right with frequency 3 GHz and unit amplitude, and the inner and outer radii of the shell 
are Ri = 0.1m and R2 = 0.2m, respectively. When Bq is equal to 0.1m and the function /(r) 
is taken as in Eq. ([8]), the radius of the virtual cylinder becomes f{Ri) = &0+-R2 = O.Sm. The 
ranges of the components of e and 71 are taken as follows: er,fir ^ [~3, — 1], eg, fig G [— 1,— |], 
and ez, fiz ^ [—3, —1]. A tiny absorptive imaginary part (~ 10^^) is added to e and JI due 
to the inevitable losses of the NRM. 

Fig. [3](a) and [3^b) are snapshots of the total electric field and scattering electric field 
induced by a PEC cylinder with radius R3 = O.Sm, respectively; Fig. [3](c) and [3](d) are 
those fields induced by the superscatterer. Comparing the patterns of electric field in the 
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FIG. 4: (a) A Gaussian beam in free space, (b) The total electric field induced by PEC cylinder 
with radius 0.6m. (c) - (d) The scattering and total fields by the superscatterer (the radius of the 
virtual cylinder is 0.6m), respectively. 



region r > R^, one can find they are almost equivalent. Here, the bounds of the amplitude 
of electric field in Fig. [3] set from —2 to 2 for clarity. The white flecks in the region of the 
complementary media show the regions where the values of fields exceed the bounds. The 
highest value of the field in the flecks is about 10^, which comes from the dominative high-m 
modes with the factor Jm{koRs)Hm {kof{r)) / Hm (koR^) in the scattering field. 

We define the magnification factor 77 as /(i?i)/i?2, which is the ratio of the radius of 
the virtual cylinder to the real size of this device. The components of permittivity and 
permeability tensors are negative for the case of r/ > 1. There is no singularity for any 
finite rj as long as both /(r) and f'{r) are nonzero in Ri < r < R2. For example, when 
Tj = 20, the ranges of the parameters are: er,^r ^ [— l-l, —0.052], eg, fig G [—19, —0.90], and 
e„/i,G [-399,-19]. 

Moreover, the device in Fig. [D^b) can be regarded as a cylindrical concave mirror for all 
angles, which can not be realized by ordinary media. From the theory of geometrical optics, 
when a plane wave incidents on it, the paraxial beams focus at r = f{Ri)/2 = {bo + i?2)/2. 



Here, a full-wave simulations with COMSOL Multipliysics are used to test the focus behavior 
of a Gaussian beam propagated from left to right with unit amplitude. Fig. 111(a)- (d) give the 
snapshots of the Gaussian beam in free space, the total electric field induced by PEC cylinder 
with radius 0.6m, and the scattering and total electric fields induced by the superscatterer 
(here -Ri = 0.1m, R2 = 0.2m and 60 = 0.4m), respectively. The scattering electric field in 
Fig. 111(c) shows the focus is approximately at r = 0.3m since Gaussian beam is close to 
a paraxial beam. Here, it is worthy noting that the superscatterer becomes a cylindrical 
concave mirror for all angles if and only if there is 69 > -^2- In principle, if we make the 
radius f{Ri) of virtual cylinder to infinity by adjusting the function /(r), the cylindrical 
concave mirror becomes a plane mirror for all angles, so that any parallel light incidents to 
the superscatterer will be reflected straight back along the incident path. 

In conclusion, we demonstrated the properties of a "superscatterer" in terms of PEC 
boundary and properly complementary media. This kind of functional devices might be 
important in EM detection. Similar concept can be extended to the case of non PEC 
boundary and three dimension. 
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